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We study the dynamical depinning following a sudden turn off of an optical lattice for a gas 
of impenetrable bosons in a tight atomic waveguide. We use a Bose-Fermi mapping to infer the 
exact quantum dynamical evolution. At long times, in the thermodynamic limit, we observe the 
approach to a non-equilibrium steady state, characterized by the absence of quasi-long-range order 
and a reduced visibility in the momentum distribution. Similar features are found in a finite-size 
system at times corresponding to half the revival time, where we find that the system approaches a 
quasi-steady state with a power-law behaviour. 
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I. INTRODUCTION 

The realization of isolated quantum systems has be¬ 
come an experimental reality with ultracold quantum 
gases. In the experiments, the atoms are confined by 
magnetic or optical means in ultra-high vacuum condi¬ 
tions and are cooled down to quantum degeneracy in vir¬ 
tually isolated conditions for a sufficiently long time with 
respect to the typical time scales for their dynamics. The 
lifetime of ultracold gases is limited to few seconds, the 
main decay mechanisms being inelastic three-body and 
spin-changing collisions. 

The study of the dynamics following a quantum quench 
in isolated quantum systems allows to address fundamen¬ 
tal questions in quantum many-body systems (see eg |T 
as well as |2 a and references therein). In one-dimensional 
integrable systems a relevant issue is the absence of ther- 
malization as oserved in the paradigmatic quantum New¬ 
ton’s cradle experiment j2j- The concept of Generalized 
Gibbs Ensemble has been introduced and developed to 
describe the state of integrable systems at long-times 

m- 

Optical lattices allow to enhance the effects of interac¬ 
tions and explore strongly correlated phases m- Bosons 
confined in optical lattices display a rich phase diagram. 
For deep lattices, ie with lattice depth considerably larger 
than the recoil energy, the bosons are described by the 
Bose-Hubbard model and display a superfluid to Mott- 
insulator transition at increasing the ratio of on-site in¬ 
teraction to tunnel energy mm- Numerical and an¬ 
alytical approaches have also been employed to study 
interaction quenches in the Bose-Hubbard model (see 
eg (iS 22). In one dimension, for an arbitrarily weak 
lattice potential commensurate with the particle num¬ 
ber density, the bosons still display a gapped, insulating 
phase, provided that interactions are sufficiently strong. 
This pinning transition is well accounted for within the 
Luttinger-liquicl description and a renormalization group 
analysis of the effect of the lattice, which yields a sine- 
Gordon Hamiltonian pTj . Various types of quantum 
quenches in the sine-Gordon model have been studied 


with different techniques eg. refermionization [22] . flow 
equations [23], and form factors [23] . 

Ultracold atomic gases confined by tight waveguides 
have been experimentally realized and the strongly in¬ 
teracting regime has been reached and characterized in 
detail j25H25] , In this one-dimensional geometry, the 
pinning transition has been experimentally observed for 
atoms subjected to a longitudinal weak optical lattice 

Ell¬ 
in this work, we consider specifically a quench across 
the pinning transition. In detail, we consider as ini¬ 
tial state a one-dimensional bosonic gas in the limit of 
infinitely strong interactions (or Tonks-Girardeau gas), 
confined by an optical lattice and with a number den¬ 
sity commensurate with the lattice spacing [3D]. We 
follow its time evolution following a sudden turn-off of 
the lattice while keeping the one-dimensional waveguide 
still present. This allows us to study the quench from 
an initially pinned, insulating ground state to an out- 
of-equilibrium depinned state. We describe the dynam¬ 
ical evolution of the system at arbitrary times using 
an exact mapping solution due to Girardeau ED E2- 
This exact solution allows us to obtain the full dy¬ 
namical solution for the quantum dynamics, going be¬ 
yond the low-energy Luttinger-liquid model or confor¬ 
mal field theory approaches. Our method allows us also 
to study the approach in time to the non-equilibrium 
steady state. Finally, we focus on the experimentally 
relavant condition of a finite-size system, choosing a ge¬ 
ometry that is amenable to experimental realization with 
ultracold atoms, eg by implementing a one-dimensional 
box-potential confinement to which an optical lattice 
is superimposed. 


II. MODEL AND EXACT SOLUTION 

We consider M bosons of mass m at zero temperature, 
confined by a longitudinal box trap of size L. They are 
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Figure 1: (Color online) Time evolution of the particle density n(L/7,t) at position x = L/7 as a function of time t 
for a commensurately filled lattice for various values of boson numbers as indicated in each panel. The horizontal 
dashed lines correspond to the (quasi) steady-state prediction in Eq. (|9j). 


described by the Hamiltonian 
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(1) 


The atoms are subjected to an optical lattice and to a 
box potential, ie V(x) = Vl cos 2 (£;lx) + V b (x), where 
kL = Mn/L and the box trap V b (x) is described by im¬ 
posing hard-walls boundary conditions on the interval 
[0,i]. The number of bosons is chosen such as to en¬ 
sure unitary filling of the lattice. The interactions among 
the atoms is modelled by contact interactions of strength 
g. In the following we shall focus on the impenetrable- 
boson or Tonks-Girardeau limit corresponding to the 
limit g — > oo. This amounts to replacing the interac¬ 
tion term in the Hamiltonian ([l]) by the cusp condition 
^(...Xj = x£...) = 0, imposing the vanishing of the wave- 
function at contact for each pair of particles {j, 1}. 

At time t = 0“ we assume that the gas is at equilibrium 
in the optical lattice. We study the time evolution of the 
gas following a sudden quench of the lattice amplitude Vl 
to zero. The exact dynamics is described by the time- 
dependent Bose-Fermi mapping, stating that the time 
evolution of the bosonic wavefunction 'F B (xi...XM>t) can 
be obtained in terms of the one of a non-interacting Fermi 
gas subjected to the same time-dependent external po¬ 
tential according to 


( 2 ) 

where ^ F (x 1 ,x 2 -,x M ,t) = det (x k , t )}. Note 
that the solution ([2]) satisfies the cusp condition at all 
times. The single-particle orbitals ipj( x k,t) are the solu¬ 
tion of the time-dependent one-body Schrodinger equa¬ 
tion 


As initial condition ipj(x, 0) we take the equilibrium 
single-particle problem in the presence of the lattice, cor¬ 
responding to a Mathieu equation with hard walls bound¬ 
ary conditions, 

H 2 

- —dlipj(x) + [V L cos 2 (k L x) + V b (x)]il)j(x) = Eji/j^x). 

(4) 

In the following it will be useful to scale all the ener¬ 
gies in units of the recoil energy Er = h 2 k\/2m , and set 
A = Vl/Er. The solution of Eq.Q is given by a gener¬ 
alization of Mathieu functions. I t am ounts to search for 
a solution of the form ipj(x) = \J2~/~LY^ n bn sin(mrx/L) 
and determine the coefficients b n . Substitution onto the 
Schrodinger equation @ yields the linear algebra prob¬ 
lem 

E sin (nk L x)(b^\n 2 /M 2 -a) + q(b^l 2M + b^l 2M )) = 0 

n 

(5) 

where a = Ej/Er — A/2, q = A/4. Equation ^ corre¬ 
sponds to an eigenvalue problem on a semi-infinite ma¬ 
trix. For the ground state of the TG gas we are inter¬ 
ested in the first M eigenvalues and eigenvectors. These 
are obtained by numerical diagonalization, performing a 
truncation to a matrix size S M. The time evolution 
of the single-particle orbitals after the sudden depinning 
is then readily given by 

ipj(x,t) = \J2/L e~ ZSrit / h sin(mrx/L) (6) 

n 

where e n = h 2 (mr/L) 2 /2m. 


III. TIME EVOLUTION OF THE DENSITY 
PROFILES 


h 2 

- -— d 2 %j}j(x, t) + V b (x)ipj(x , t ) = ihd t ipj(x, t), (3) 


In order to explore the post-quench dynamics we ana¬ 
lyze the time evolution of various observables. We con- 
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sider first the time evolution of the density profile. This 
is obtained, using the Bose-Fermi mapping, as the one of 
the corresponding Fermi gas, 

M 

n(x,t) = (?) 

3 

Substitution of the explicit solution for the lattice prob¬ 
lem yields 

M 

n(x,t) % 

j n,n' 

xsm(mrx/L)sm(n'nx/L). (8) 

The time evolution of the density is shown in Fig. [T] for 
various numbers of bosons and various fixed sizes, at con¬ 
stant filling of one boson per site. Recurrences are clearly 
visible, as expected for a finite-size system. While one 
could estimate as trivial recurrence time T = 2itH/ei, we 
note that the density profiles show revivals at earlier time 
Tf> = T/AM. This property is specific to our choice of 
system and initial state. The time evolution in Eq. © is 
determined by the energy difference e n —£ n ' oc (n 2 — r? 2 ); 
the most important contribution in the coefficients is 

for n' — n = 0 mod 2M, as one can infer from the nu¬ 

merical solution of Eq. © as well as from a perturbative 
approach J3D] at weak lattice strength, thereby yielding 
the observed recurrence time Tr. 

IV. THE NON-EQUILIBRIUM STEADY STATE 

For a value of M sufficiently large to be outside the 
mesoscopic regime of a few particles (M > 11 in our case) 
and for sufficiently long times (ie at about half revival 
time) we observe in Fig. [l] that the system tends to a 
(quasi) steady state - as we shall denote this state for a 
finite-size system. |34| This state tends to a truly steady 
state once the thermodynamic limit is taken (see Section 
[Vjbelow for details). The (quasi) steady state can be well 
described by neglecting the oscillating terms in Eq. © 

2 °° 

n(x, t) —> n ss (x) = — f n sin 2 (nnx/L) (9) 

1j ^ J 

n—1 

with nonthermal occupation numbers /„ given by f n = 
HjLi \bn^\ 2 - The above result M is illustrated as hor¬ 
izontal dashed lines in Fig jT] In “the absence of the lat¬ 
tice f n is given by a Fermi distribution at zero tempera¬ 
ture. At increasing height of the initial lattice the (quasi) 
steady state is characterized by the occupation of more 
and more excited bands, as shown in Fig. [2] We notice 
that the occupation numbers f n vanishes for n = M + 1 
to 2 M, and similarly for higher excited levels correspond¬ 
ing to even bands. This is due to the fact that the op¬ 
tical lattice acts as a backscattering potential creating 
excitations with wavevector 2fci, which corresponds to 
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Figure 2: (Color online) Occupation numbers f n as a 
function of the quantum number n for M = 50 
particles, for various values of the dimensionless lattice 
strength, from left to right, A = 0,4, 50. 

2kp = 2 Mtt/L in our choice of lattice filling. In en¬ 
ergy space, this allows to excite only levels with quantum 
number difference An = 2 M. Mathematically, the result 
follows from the linear algebra problem in equation (5) : 
the optical lattice gives rise to off-diagonal terms into the 
matrix that are a distance 2M from the diagonal. The 
corresponding eigenvectors have mostly zero components 
except bn 1 0 ,Vn = 2 M + 1. Combining this prop¬ 
erty with the definition of the occupation numbers f n we 
obtain the result shown in Fig.2. 

Furthermore, our exact solution allows also to explore 
the approach to steady state. By an extensive analysis 
of various systems sizes, we have found a power-law ap¬ 
proach to the steady state ie, n(x,t ) — n ss {x ) ~ 1/t as 
times approaches Tr/ 2, which corresponds to the large¬ 
time limit in our finite-size system. This is illustrated in 
Fig|3]for various values of particle number M. 

V. FIRST-ORDER COHERENCE AND 
MOMENTUM DISTRIBUTION 

In order to further characterize the properties 
of the (quasi) steady state we study the time- 
dependent one-body density matrix pi(x,y,t) = 
M f dx 2 --dx N 'l’* B (x,X 2 , ...x M ,t)'f> B (y,x 2 , This 

allows to determine the coherence properties of this state 
and in particular the presence of quasi-off-diagonal long- 
range order (QODLRO). Furthermore, this allows to ob¬ 
tain the momentum distribution of the gas n(fc, t) = 
f dx f dypi(x,y,t)e lk ^ x ~ v \ which is experimentally ac¬ 
cessible with a high precision (see eg [55]). 

Following the approach of Ref. ESI the one-body den¬ 
sity matrix of a time-evolving Tonks-Girardeau gas is 
given by 

M 

Pi(x,y,t)= ^ -i/i*(x,t)Aji(x,y,t)ipi(y,t). (10) 

3 , 1 =i 
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Figure 3: (Color online) Approach to nonequilibrium 
steady-state: time evolution (time in units of Tr) in 
double logarithmic scale for the difference between the 
particle density and its corresponding steady-state value 
\n(x,t) — n s,s (a:)| (in units of no = M/A), evaluated at 
x = A/7 and for A = 40, for 81 (blue), 111 (red), 151 
(orange) and 201 particles (purple). The particle 
density has been time averaged over a short time 
interval At = 0.001 Tr in order to decrease the noise in 
the figure. The dashed line indicates the 1/t power law 
decay. 


where the matrix A(x,y,t) = (P 1 ) T detP, 

P(x,y,t ) = 1 - <5, with Qji(x,y) = 2sign (y - 
x ) fj'dx'i/>j(x',t)i/>i(x',t)- For our specific case, the ma¬ 
trix elements Qji(x,y) are readily evaluated analytically 
using Eq. |6|. 

The resulting density matrices are illustrated in Figj4] 
for 15 bosons. As compared with the equilibrium case 
in absence of the lattice (panels a) and d)), the effect 
of the lattice at time t = 0 is a pinning along the di¬ 
agonal x = y and and a reduction of the off-diagonal 
coherences, as evident from the sections taken along the 
direction x = —y (panels e) and f)). As a main result, 
we find that the (quasi) steady state at time t = Tr/2 
displays no QODLRO, and the one-body density matrix 
decays exponentially EZ1 This is especially striking since 
the system is evolving in a homogeneous box, ie it is de- 
pinned: While the corresponding equilibrium state in the 
box displays QODLRO, ie the well-known power-law de¬ 
cay with power exponent -1/2, as shown in panel d) of 
Fig. a quite remarkably, the non-equilibrium aspect of 
the gas influences dramatically its coherence properties. 
This feature is found quite generally in integrable mod¬ 
els, and has been first predicted by conformal-field meth¬ 
ods for the state of the system at long times following a 
quench across a quantum critical point [I]. 

For weak lattice strength, the exponential decay of the 
one-body density matrix for the (quasi) steady state can 
be analytically obtained: the main contribution to the 


weights 5/, is given by the term S n j , yielding for the ma¬ 
trix elements Pij = dij[l — 2\(x — y)/L—[sm(2wjx/L) — 
sin(27rj2//A)]/(27rj)|]. At large relative distances one may 
set Pij ~ Sij( 1 — 2\x — y\/L), and thereby obtaining 
A = 1(1 — 2\x — y\/L)( N b. Taking the thermodynamic 
limit M — > oo and A —► oo at fixed n = M/A we find 
A = 1 exp(—2n|a; — y |) and finally using (101 we obtain 
the bosonic quasi-steady state one-body density matrix 

pi{x,y,t~> oo) = p 1F (x,y,t —> oo) e~ 2n \ x - y \ , (11) 


where p 1F (x,y,t) = t) is the 

fermionic one-body density matrix. This result coin¬ 
cides with the one obtained in m for the release of 
a harmonically trapped TG gas onto a ring: as ex¬ 
pected, boundary conditions do not affect the result in 
the thermodynamic limit. Notice however the different 
way the steady state is reached: while in Ref. [T3] the 
long-time state is obtained by summing up all the pe¬ 
riodically repeated images on the ring, in our case it is 
obtained by multiple reflections at the boundaries. No¬ 
tice that our result can then be linked to the concept 
of Generalized Gibbs ensemble (GGE) for the thermo¬ 
dynamic limit of our model: as in M, the average of 
the fermionic occupation numbers (c^c*,), obtained from 
the initial fermionic one-body density matrix according 
to (cj)cfc) = f dxf dyp 1F (x, y, 0)e lk( - x ~ v ' ) /(2 n) are con¬ 
served in the time evolution and can be used to determine 
the Lagrange multipliers A & defining the density matrix 
of the system at long times: pgge ~ exp(— )T) fc A^c^c*,). 

From the knowledge of the bosonic one-body density 
matrix pi(x,y,t ) we finally obtain the exact momentum 
distribution of the gas. We stress that our exact approach 
allows us to cover all the ranges of momenta, beyond 
the low-momentum region accessible by conformal-field 
methods. As shown in Fig(5j the momentum distribu¬ 
tion of the (quasi) steady state is considerably different 
from the equilibrium ones both in absence and in pres¬ 
ence of the lattice: it displays a considerably reduced 
intensity at low momenta, and does not show the typical 
backscattering peak at k = 2k F found for the equilibrium 
gas in the presence of the lattice. 


VI. CONCLUSIONS AND OUTLOOK 

In conclusion, we have studied the exact time evolution 
of a Tonks-Girardeau gas following its sudden depinning 
off a weak optical lattice. We have identified a suitable 
long-time limit where a non-equilibrium steady state is 
reached in the thermodynamic limit, and we have shown 
a power-law approach to the steady state. Furthermore, 
we have shown that this state is characterized by the ab¬ 
sence of quasi-long-range order ie an exponential decay 
of one-body correlations, in agreement with the predic¬ 
tions of the Generalized Gibbs Ensemble. Our numer¬ 
ical analysis for a system of finite size shows that this 
scenario could be reached with experimentally realistic 
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Figure 4: (Color online) Real part of the one-body density matrix pi(x,y) in units of no = M/L as a function of the 
coordinates of the spatial coordinates x and y for M = 15 bosons. Top panels: top view, bottom panels: 
corresponding cuts at fixed x + y = L /2 as a function of the relative distance \x — y\. a) and d) equilibrium state in 
absence of the lattice, b) and e) equilibrium state in presence of the lattice, with lattice strength A = 10. c) and f) 
depinned quasi steady state at time t = Tr/2 after a quench of the lattice to zero. The dashed line in panel e) is a 
guide to the eye and in panel f) indicates the exponential decay e ~ 2n ' x ~ v I predicted in Eq.(11). 
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Figure 5: (Color online) Momentum distribution n(k) 
as a function of wave vector k for M = 15 particles. 
Left: equilibrium state in absence of the lattice, center: 
equilibrium state in presence of a lattice with A = 10, 
right: quasi steady state at time t = Tr/2 after a 
quench of the lattice to zero. 


numbers of bosons in a tight atomic waveguide, and that 
the time-dependent momentum distribution yields rele¬ 
vant information about this state. Our work opens to the 
study of the details of the quench depinning dynamics of 
a Lieb-Liniger gas at arbitrary interactions and to further 
tests of the GGE hypothesis with hard walls boundary 
conditions (see eg |3B]). 
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